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The initial mass spectrum of old globular clusters in dwarf galaxies 
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ABSTRACT 

We test whether the masses of old globular clusters (GCs) in dwarf galaxies are consistent 
with the same initial mass spectrum as young massive clusters (YMCs) in nearby star-forming 
galaxies. The most massive GCs of dwarf galaxies are compared to their expected masses 
when drawing from the Schechter-type ICMF of YMCs. It is found that the most massive 
GCs of galaxies in the stellar mass range M^ gai = lO^-lO'' Mq are consistent with the same 
initial mass spectrum as YMCs in about 90% of the cases, suggesting that their formation 
mechanisms were the same. For the remaining 10%, the most massive clusters are nuclear 
GCs, which have been able to grow to higher masses through further merging after their ini- 
tial formation ended. Because the effects of cluster disruption are weaker for more massive 
clusters, we estimate that up to one third of the metal-poor GCs in the Milky Way may have a 
nuclear origin, while the remaining two thirds formed through the same process as YMCs in 
the local universe. A log-normal ICMF is inconsistent with observed GCs at a 99.6% confi- 
dence level. 

Key words: galaxies: evolution - galaxies: dwarf galaxies - galaxies: starburst - galaxies: 
star clusters - globular clusters: general 



1 INTRODUCTION 

The quest to understand how globular clusters (GCs) formed goe s 
back several decades (e.g. lSearle & Zinij|l978l : lFall & Reej|l985h . 
Recent efforts have concentrated on the formation sites and mech- 
anisms of GCs, driven by the question whether GCs are presently 
still forming in the nearby universe. The Gal actic GC population 
. likely has its roots in several e nvironments jMackev & Gilmor3 
■ |2004 iMuratov & GnedinllZOlOh . Metal-poor GCs are thought to 
, originate from dwarf galaxies that were accreted by the Milky Wa y 
jPrieto & G nedin 2008; Pefiarrubia et al. 20oA iLee et"al] |2009|) . 
while metal-rich GCs may have formed in gas-rich galaxy merg- 



young cluster populations. Indeed, in nearby spiral galaxies, star- 
bursts and galaxy mergers, the masses of young massive clusters 
(YMCs) ar e all consistent with a n approximately univ ersal ICMF 
dCieles etal. 2006; Larsen 2009i; IPortegies Zwart eTal. 2010). ex- 
cept for a n environmentally d ependent variation of its upper 
truncation jSastian et al.l l20lTbh . The approximately log-normal 
present-day mass function of GCs also bea rs some traits of forma- 
tion according to a Schechter -type ICMF i Harris & Pudritz 1994|; 
McLaughlin & Pudritz 199^; iBurkert & Smi th 200Q; Jordan et all 
20071), but this necessarily implies that the vast maj ority of low- 



ay ■ 

ers ^Ashman & ZeDjl992l;IWhitmore et all 1 9991 ; IZepf et al J 19991) 
or unstable galaxy discs i Shapiro et al.ll20ld ). Whichever environ- 
ment GCs have formed in, their masses suggest it should have been 
in a burs t of star formation, potentially similar to ne arby galaxy 
mergers faoltzman et al. 19921; Schweizer et aljn"996l) or starburst 
dwarf galaxies jAdamo et alj|2010l) . 

A starburst origin for GCs is required if their initial clus- 
ter mass function (ICMF) was the same as that of young star 
cluster population s in th e local universe, which follows a power 
law or ISchechteil j 19761) t ype distribution with index —2 down 
to some lower mass limit jPortegies Zwart et al.ll201(il) . The ran- 
dom sampling from such an ICMF requires a high star forma- 
tion rate (SFR) for the production of massive clusters, unless the 
lower mass limit exceeds ^ 10^ M©, which is not consistent with 
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mass GCs has been destroye d due to tidal disruption /Fall & Reel 
1977; Elmegreen & Efreinq^ 



Il997l; |\^sperini 2001; Fall & Zhan 
20011 ; IKruiissen & Portegi es Zwartll2009l) . This scenario contrasts 



with the idea that the current mass function of G Cs would reflect 
a fundamentally different (e.g. log -normal) ICMF iVesperini et al.l 
l2003l ; |Parmentier & Gilmor3 2007), which has been put forward as 



an alternative explanation for the approximate universality of the 
peaked present-day GC mass function. A siinilar ICMF for GCs 
and young clusters in nearby galaxies would imply that the dis- 
ruption histories of GCs have also been close to universal. Consid- 
ering the broad range of galactic environments in which GCs are 
currently found, this suggests that most of the GC disruption may 
have occurred at the epoch of their formation dElmegreenl |20 1 ol ; 
IKruiissen et alj20"Tll Kruijssen et al. in prep.), when their environ- 
ments were likely more comparable. 

It is relevant to ask whether or not the ICMFs of GCs and 
young clusters in the local universe are the same, because the ICMF 
is a tr acer of the cluster formation process jElmegreen & Falgarond 
Il996h . A similar ICMF would thus be required if GCs and YMCs 
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are two stages of the same type of object in terms of their physics 
and formation mechanism. Such evidence cannot be found for the 
entire GC mass range, due to the ongoing disruption of (low-mass) 
stellar clusters, but this problem can be addressed using the most 
massive GCs. Given an ICMF, the mass of the most massive star 
cluster can be predicted statistically if a fraction of star formation 
resulting in bound clusters is assumed, i.e. the cl uster formatio n 
efficiency or CFE. This has been done for YMCs teastian 2008'), 
but for GCs such an approach has been obstructed by unknowns 
about the conversion of current GC properties to their original form. 

We aim to connect YMCs and GCs by adopting the recent ad- 
vancements made in the understanding of YMC formation in the 
local universe. These will be used to predict the expected masses 
of the most massive GCs in dwarf galaxies for limiting cases, to 
verify whether they are consistent with random sampling from 
the same ICMF as YMCs. This should provide an indication of 
whether or not their formation mechanisms are compatible. Dwarf 
galaxies are the most suitable for such an analysis, because they 
have not accreted as much since GC formation as massive galax- 
ies have. Moreove r, the masses of the most massive YMCs (e.g. 
ISastian et alj |2006') exceed the stellar masses of the faintest dwarf 
galaxies (e.g. (Mateo 1998), and therefore dwarf galaxies are the 
prime targets in which the limits of GC formation can be explored. 
This naturally constrains our analysis to the origin of metal-poor 
GCs. 



2 THE RELATION BETWEEN GLOBULAR CLUSTER 
MASS AND GALAXY STELLAR MASS 

Young star clusters in the Milky Way and nearby gala xies follow 
a pow er law ICMF with an exponential truncation (cf. Schechtei' 



\91t ) , over a mass rang e of M = 10-10" Mfl (M ai-aston e t al. 
2OO4I : iGieles et al.ll2006l ; IPortegies Zwart et al.llS)Iol) . The proba- 
bility density of cluster masses ^(M) is expressed as 



^(M)M oc M"^ exp(-M/Mc)M, 



(1) 



where Mc represents the truncation mass and ^(M) is normalised 
such that J ^(M)dM = 1, integrated between the physical mass 
limits Mmin and Mmax. In this study, we adopt Mmin = 10^ Mq. 
The physical upper mass limit is set equal to the stellar mass of the 
host galaxy, i.e. Mmax = M*,gai, with generally Mc < Mmax- 

When randomly sampling cluster masses from an ICMF, the 
probability distribution function for the most massive sampled clus- 
ter p(M„,ax) is given by jMaschberger & Clark3l2008h : 



p(M„ax) = JV 



C(M')dM' 



C(M„,ax), 



(2) 



where N is the total number of clusters. It is obtained by integra- 
tion of Eq.[T]after renormalising the ICMF to match the total stellar 
mass formed in clusters M*_ci. The expected Mmax is given by the 
median of Eq|2l with the 1 6th and 84th percentiles spanning one 
standard deviation up and down. We use Eq.|2]to compute the most 
massive cluster as a function of host galaxy stellar mass, and com- 
pare the result to observed GCs and (for verification) YMCs. 

The comparison between the theoretically expected most mas- 
sive GC and observed GCs requires certain assumptions about the 
following quantities and their dependences on the host galaxy. 

(I) The fraction of the present-day galaxy stellar mass that was 
formed coevally with the GCs (f*,oid). 
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Figure 1. Mass of the most massive GC versus the stellar mass of its 
host galaxy. The black solid line shows the theoretical upper limit (see 
text) to the relation for a power law ICMF. with the standard devia- 
tion due to random sampling from Eq. |2] indicated as dashed lines. Grey 
lines denote the same for a log-normal ICMF (see text). Coloured sym- 
bols denote observed GCs, with host galaxy types indica ted by the leg- 
end. Squares mark GCs from nearby field dwarf galaxies (Georgiev et^lJ 
2010), dot s denote GCs from galaxies i n the Virgo Cluster ( Cote et alj 
I2OO4 iPeng et alj|2008t Ijordan et al. 2009), and diamonds indicate (in or- 
der of incr easing M^. gal) Forna x I Mackey & Gilmore 2003c ; Lokas 2009^ 
the SM C iMackev & Gilmore 2003a) and the LMC iMackev & Gilmo3 
l2003bl) . Black circles mark nuclear GCs, defined as having a projected 
galactocentric radius < 0.5 kpc. The typical error margin on the data is 
shown in the top left corner. 



(2) The fraction of star formation that yields bound stellar clus- 
ters, i.e. the CFE (f). 

(3) The exponential truncation mass of the ICMF (Mc). 

(4) The mass loss suffered by GCs since their formation, 
parametrized by the disruption time-scale of a lO' M© cluster (ts). 

The parameters ft,oid and V relate the current stellar mass of the 
galaxy M*,gai to the normalisation of the ICMF through M,,ci = 
rf*.oidM,,gai, and determines what the present-day mass of the 
GCs is, after a Hubble time of dynamical evolution (see below). 

The hypothesis that old GCs and YMCs form from an ICMF 
with a common functional form can be tested to first order by set- 
ting the above quantities to values that enable the formation of the 
most massive clusters. Such a limiting case would assume that (I) 
all present stellar mass of the galaxy was formed at the same time 
as the GCs (f,,oid = 1), (2) that all of its stars formed in clusters 
(r = 1), (3) that the ICMF is a pure power law with an upper 
mass limit equal to the stellar mass of the galaxy (Mc = 00 and 
Mmax = M*„ai), and (4) that there has been no mass loss by dy- 
namical processes (ts = 00). Stellar evolution would still apply and 
is included as a fixed, 30% decrease of the cluster mass over a Hub- 
ble time. The corresponding parameter set is shown as the 'extreme 
mode' in Table [T] and gives an upper limit to the relation between 
the most massive GC and the stellar mass of the host galaxy, pro- 
vided that the GCs are randomly sampled from a lSchechten ( 1 1 9761) 
type ICMF. If any of the observed GCs lies above this relation, it 
would imply that at least some GCs are too massive to have been 
sampled from the ICMF of YMCs, and would suggest that they 
formed through a different mechanism. 

The 'extreme mode' relation between the most massive GC 
and host galaxy stellar mass is shown in Fig. [T] together with 
the observed most massive GCs from nearby field dwarf galaxies 
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and galaxies in the Virgo Cluster. For reference, we have also in- 
cluded the relation for a purely log-normal IC MF, with a peak mas s 
log(Mpeak/M0) ~ 5.1, dispersion a ~ 0.6 ( Vesperini et al.'2003'), 
and a weak dependence of both on galaxy mass as in Jordan e t al,, 
( I2OO7I) . Figure[T]shows that even in this limiting case a log-normal 
ICMF is not consistent with the observed GC masses. This is 
mainly caused by the steep decline of a log-normal function at the 
high-mass end. By contrast, none of the observed GCs has a mass 
higher than the theoretical upper limit for randomly sampling the 
most massive GC from a Schechter-type ICMF. This shows that the 
initial, total star-forming potential of the galaxy could indeed 'sup- 
port' the formation of such massive GCs through the initial mass 
spectrum we see in the local universe. It does not yet imply that the 
ICMFs of GCs and YMCs are fully consistent, because this would 
assume that the adopted parameter set represents the conditions in 
real galaxies. 

However, the scenario sketched as the 'extreme mode' of GC 
formation in Fig.[T]is not quite realistic. Only a certain fraction of 
the stellar content of dwarf galaxies formed coevally with GCs (e.g. 
iMcOuinn et al .l201(]h . there is no evidence for such efficient cluster 
formatio n that V = 1 (e.g.'God dard et al.l2OI0l : lBressert et alJlOld : 
[Silya -V illa & Larsen 201 1), there is likely an enviroi mientally de- 
pendent upper mass limit to cluster formation (e.g. iGieles et al.l 
l200dlBastiai]|2008h . and GCs are known to have lost at least 
some fraction of th eir initial mass due to dynamical evolution 
( lFall&Zhan"3 l200l'). The values of these parameters and their re- 
lation to the properties of the host galaxy have been constrained 
empirically for cluster formation in the local universe. We adopt 
these relations to test whether GC formation has proceeded in a 
way that is consistent with the ICMF of young clusters in nearby 
galaxies. 

In a st u dy o f seven nearby star-forming galaxies, 
iGoddard et all bOld) related V to the star formation rate 
density Tsfr: 



r = 0.29 



(3) 



^M© yr~' kpc" 

which has been confirmed bv lAdamo et al ]| |20ld. l20T l|). The rela- 
tion holds over E sfr = 6 x 10~^ -1 M© yr~' kpc~^, corresponding 
to r = 0.08-0.29. iBastianI l l2008i) explored the relation between the 
brightest cluster in a galaxy and the star formation rate. If the CFE 
is only allowed to vary slowly as in Eq. [3] then Mc must increase 
with the SFR. Based on his results on YMCs, we assume 



Mc = 7 X 10' 



SFR 



Mr: 



(4) 



M© yr- 

which reproduces the typical Mc ~ 2 x lO'' M© found bv iLarseiil 
( I2OO9) for spiral galaxies with SFR ~ 0.1 M© yr~' and Mc ~ 3 x 
10^ M© for the An tennae galaxies, assuming SFR = 20 M© yr~' 
dZhang et alj200ll ). Equations[3]and|4]depend on the galaxy stellar 
mass M«,gai through SFR(M*_gai) and IIsFR(M*.gai) (see below). 

After their formation, star clusters lose mass due to stellar evo- 
lution and tidal disruption. While this may not be very important for 
YMCs, the mass of GCs will have decreased after almost a Hubble 
time of evolution. A simple empirical relation for the mass evolu- 
tion of a cluster was derived by Lamers et al. ( 2005a) : 



M(t) = M, 



/isc 



Mi 



105 



1/7 



(5) 



which depends on the initial c luster mass M , the s tellar evolution- 
ary mass loss (^sc ^ 0.7, e.g. lLamers et al.l 2005ah . the parameter 



7, which sets the mass -dependence of the disruption time (we adopt 
a value of 7 ~ 0.7, see lLamers et alj20IOl) , and the disruption time- 
scale of_a_J^^_^l^cluster tj, which is environmentally dependent 
(e.g. ' Kruiissen et al.ll20I ll : iBastian et al. 201 la). In their analysis 
of GCs in the Virgo Cluster, Jordan et al] j2007l) find that the peak 
of the GC luminosity function (GCLF) depends on the magnitude 
of the host galaxy. By assuming that the peak lum inosity is a pprox- 
imately proportional to the disruption time-scale tGieles 2009) and 
by using a constant M/L ratio to convert the galaxy magnitude to 
its total stellar mass, this can be translated to a very weak depen- 
dence of ts on Mt.gai : 



Cdi 



M* 



gal 



10' M, 



(6) 



where log(Cdis/yr ) ^ 10.5 for SMC /LMC-type galaxies with 
M*.gai - 10' M© jLamers et aLll2005bh . For the Milky Way, this 
gives t; ~ 2 X 10'° yr, consistent with the range derived for Galactic 
GCs iKruiissen & Mi eske 2009). The scatter among GCs in a par- 
ticular galaxy is non-negligible though, and we therefore assume a 
spread of a factor of two, i.e. log (Cjis/yr) = 10.2-10.8. 

By combining the above set of equations with Eq. |2l we can 
relate the galaxy stellar mass to the most massive GC that is ex- 
pected if the GCs were formed with the same ICMF as YMCs in 
the nearby universe. For this, we adopt certain relations or values 
for SFR(M*,gai), 5IsFR(M*,gai), f*,oid and Cdis. Due to the depen- 
dence on the SFR and Esfr, we need to define two limiting modes 
of GC formation that span the range of plausible star formation 
histories ( SFHs) of the gala xies. We use the semi-analytic galaxy 
mode ls of FGuo et all lioTH). based on the Millennium-II simula- 
tiorQ teovlan-Kolchin et al.l2009D . This then spans the range in the 
Mmax-M*,gai plane that would be expected if GCs were formed ac- 
cording to the ICMF of YMCs. The two limiting modes are: 

(1) The 'quiescent mode', which assumes that the SFR has been 
constant throughout the history of the galaxy. In this mode, f*.oid is 
set by the definition that any star or cluster formation corresponding 
to ages r > 5 Gyr is considered to be coeval with GC formation. 
While this age limit is arbitrary to some degree, it ensures that we 
are considering the galaxies at times when GCs may ha ve formed. 
We ha ve fitted IIsFR(M*,gai) by selecting galaxies from lGuo et al.l 
( l201lb at the moments when they are forming stars at a rate lower 
than their mean SFRs. This selection guarantees that any episodes 
of enhanced star formation are excluded (as, by definition, the SFR 
during these episodes will be greater than the average SFR over 
the history of the galaxy). For a given galaxy, Esfr is obtained by 
dividing its SFR by the area within its half-mass radius at each 
quiescent time-step, and then taking the average of these values. 
This approach accounts for the relation between the SFR and Esfr 
as the structure of the galaxy changes. 5Isfr(M* „ai) is then obtained 
by fitting a power law relation to the galaxy sample. We also adjust 
Cdis to the low end of its range of possible values. This mode leads 
to the lowest maximum GC masses. 

(2) The 'rapid burst mode', which assumes that all stars in the 
galaxy were formed coevally with the GCs in a starburst. This 
implies M*_gai = Mburst, with Mburst the total mass formed in the 
burst. We have fitted S FR(Mburst) by a gain selecting galaxies at 
ages r > 5 Gyr from IGuo et al] ( 1201 ll) . this time picking star- 
bursts (i.e. the times at which the SFR of each galaxy peaks) with 
durations shorter than the time resolution (300 Myr). We then as- 
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Table 1. Adopted parameters. For any SFR and 5Isfr> Eqs.[3]and|4]give the values of V and Mc. Only their resulting dependences on M* gai are listed here. 



SF mode 


^*,old 


SFR (Mq yr-') 


5:sFR (Mq yr ' kpc ^) 


r 




Cdis (yr) 


Extreme 


1 






1 


oo 


oo 


Rapid burst 


1 


32.9(M*.gi,l/10'' Mq^-"^^ 


4.53(M^_gal/10'' Mq)"-* 


(X m0-'4, 
*.sal 




1010.8 


Mild burst 


0.8 


0.28(M*,gai/10' Mo)"-'*" 


0.033(M*.gal/10'' Mq)"-*' 


*,gal 


oc M"'^ , 


1010.5 


Quiescent 


0.6 




0.021(M*,gai/10'' Mq)'"5 


*,gal 


OC M°-^ , 


1010.2 



sume burst durations of one dynamical time to determine the actual 
SFR, i.e. tburst = 3R/V with R the disc scale length and V the ro- 
tational velocity. Again, the area covered by the half-mass radius is 
used to determine Hsfr- We adjust Cdis to the high end of its range 
of possible values. This mode leads to the highest maximum GC 
masses. 

The resulting parameters and proportionalities for r(M*.gai) and 
Mc(Mj,_gai) are shown for all scenarios in Table [T] In addition to 
the above two modes of GC formation, we include an interme- 
diate, 'mild burst' scenario for which we have fit ted SFR(M,t,gai) 
and 5IsFR(M*_gai) to galaxies from lGuo et al. at the moments 

when they are forming stars at a rate higher than their mean SFRs, 
analogously to the quiescent mode above. This selection ensures 
that only episodes of enhanced star formation are included. 



3 THE FORMATION OF GLOBULAR CLUSTERS 

The region in the M*_gai~Mmax plane that can be spanned by random 
sampling from a Schechter-like CIMF for the physically plausible 
range of dwarf galaxy histories is shown in Fig. [2] The top panel 
shows the relation at the time of GC formation and includes the 
most massive YMCs of nearby galaxie^ for reference. For these 
data points, the position on the x-axis does not reflect the total stel- 
lar mass of the galaxies, but merely the mass that was formed co- 
evally with the most massive cluster - either given by the integral 
of the SFR over the starburst or, if the SFR has been constant or 
the SFH is unknown, the product of the current SFR with twice 
the age of the YMC. The YMCs are scattered along the top edge 
of the 'allowed' region, because most of them reside in starburst 
galaxies and we have implicitly assumed f*_oid = 1 by using the co- 
eval stellar mass. Acknowledging this, we see that the YMCs fall 
in the anticipated range for the quiescent and rapid burst modes, 
which illustrates that our approach is indeed consistent with their 
formation. Due to their Schechter-type ICMF, the YMCs are not 
consistent with the relation for a log-normal ICMF. 

The bottom panel of Fig. [2] shows the observed and predicted 
present-day masses of the most massive GCs in the dwarf galaxies 
of our sample as a function of the total host galaxy stellar mass. In 



In order of increasing coeval stellar mass, the galaxies denoted by 
the symbols are NGC1705 (ref. 1,2), Mrk930 (ref 3), M83 (ref 4,5), 
NGC4214 (ref 5,6), NGC1569 (ref 1,7), NGC6946 (ref 4,5), NGC4449 
(ref 5,6), NGC4038/39 (ref 8,9,10), NGC3256 (ref 11,12), Haroll (ref 
13,14) and M82 (ref. 15,16). The sources that were used to c ompil e 
the data are: JJ) I Smith & Gallagh ej foOOll), (2 ) lAnnibaU et alj j2003| 



(3) lAdamo eTai] 4201 ih, (4) iLars en & Richtl eil )2000l). (5) iLarsen et alj 
i2004h. (6) iMcOuinn et"ai] 1201 q). ( 7 ) .Grocholski et alj j2008l), (8 ) 
IZhang etalJ )200lh. (9 ) iMengel et ai] )2002h, (10) iKarl et alj j2010l) 
(11) IZepfetal.' ('l999|), (12) j Trancho et al.' fe00 7|), (13) iHaves et alj 
(2007), (14) Adamo e t all J2010l) ~(15) Forster Schreiber et alj i l2003l) . (Ti) 
iKonstantopoulos et alj i2009l) . 
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Figure 2. Mass of the most massive star cluster versus the stellar mass of 
its host galaxy. Solid and dashed lines indicate the upper Hmit and its stan- 
dard deviation from Fig.[T] The shaded area represents the parameter space 
encompassed by the 'quiescent mode' and 'rapid burst mode' from Table[T] 
which represent the physically plausible extremes. The region includes one 
standard deviation up and down. The long-dashed line indicates the 'mild 
burst mode'. The dotted lines mark the equivalent of the shaded area for 
a log-normal ICMF (see Sect.|2j. Top: initial relation, excluding mass loss 
due to stellar evolution and tidal disruption. Symbols indicate YMCs in the 
nearby universe, coloured by their ages r. Their location on the x-axis is 
determined by the total stellar mass formed coevally with the YMC (see 
text). Bottom: relation at the present day, after the application of mass loss 
due to stellar evolution and tidal disruption. Symbols represent the sample 
of GCs and their host galaxies as in Fig.[T] The typical error margins on the 
data are shown in the top left comers of each panel. GCs that are more than 
2cr separated from the 'allowed' area are marked with red boxes, which are 
connected by red lines to the second most massive GCs in these galaxies. 



this panel, mass loss due to stellar evolution and tidal disruption has 
been included, which skews the region that is spanned by the quies- 
cent and rapid burst modes and moves it to lower masses. There is 
good agreement with the observed GCs, suggesting that the forma- 
tion of GCs is indeed consistent with the same ICMF as YMCs. It 
also indicates that the CFE and Mc followed the same relations for 
GCs as they do for YMCs. If the ICMF were not (exponentially) 
truncated, most of the Virgo Cluster galaxies would fall well below 
the predicted region. 
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At least some discrepant GCs are expected statistically, and 
indeed a handful of GCs have masses that are Icr to 2(j outliers 
from the predicted range. For the two dSph galaxies at M*,gai ~ 
2 X 10'' Mq, this could have a physical explanation, such as a very 
efficient starburst that formed the entire galaxy at a CFE of f = 1. 
Looking at statistically more substantial extremes, we see that three 
out of the 101 GCs have masses that are more than 2a separated 
from the 'allowed' region in the parameter space (marked by red 
boxes in Fig.|2]l, whereas a purely log-normal ICMF is more than 
2(T-inconsistent with the masses of 16 GCs, and generally cannot 
reproduce the observed spread in the Mmax-Mt.gai relation. When 
leaving out nuclear GCs (nGCs, see below), a total of eight 2(7- 
inconsistent GCs are left, ruling out a log-normal ICMF for metal- 
poor GCs in dwarf galaxies at a 99.6% confidence level. By con- 
trast, there are no non-nuclear outliers at all for the Schechter-type 
ICMF. We have tested the agreement for log-normal ICMFs with 
different parameters and find that a log-normal can only yield suffi- 
ciently hi gh GC masses for disp ersions cr > 1 .4 (which is inconsis- 
tent with Vesperini et al.ll2003f) , but this would still not reproduce 



the large spread in Mmax. The possibility of a higher peak mass 
is ruled out by current GC systems, because dynamical evolution 
is not capable of then reducing t he peak mass to its current valu^ 
jjordan et alj|2007l : lGielesll2009l) . 

All three GCs that are more than 2(T-inconsistent with sam- 
pling from a Schechter-type ICMF are nGCs (black circles in 
Fig.|2]l. In these galaxies, the second most massive GCs are con- 
sistent with the predicted mass range as well as with the masses of 
the GCs in the other galaxies. Their masses require them to have 
formed in the rapid burst mode, especially since they are only the 
second most massive GCs. The abundance of nGCs among the 
outliers is not necessarily surprising, since nGCs are thought to 
be the products of a combination of merged star cl usters that mi- 
grated to the galaxy centre due to d ynamical friction jMiocchi et al.l 
I2OO6I) and additional gas accretion iHartmann et al.ll201 1 ). This al- 
lows their mass to increase beyond the limit imposed by the fi- 
nite duration of the star formation process. The form of the ICMF 
is the outcome of the hierarchical m erging of star-forming aggre- 
gates telmegreen & Falgarone|[l996l) , and is therefore normalised 
to the total mass forming in stellar clusters. By contrast, the nGCs 
are consistent with the 'extreme mode' of random sampling from 
Fig. [T] in which the stellar mass of the entire galaxy is converted 
into an ICMF. They could thus be regarded as the high-mass end 
of the mass spectrum of all (stellar) structure in the galaxy, driven 
by hierarchical merging over a Hubble time rather than on a star 
formation time-scale. 

The fraction of galaxies with a 2(T-outlier as their most mas- 
sive GC is ~ 10% between M*,gai = lO'-IO' Mq , and decreases for 
more massive ga laxies. Since nGCs ar e likely formed through dy- 
namical friction teellazzini et alj2008l) . this is probably caused by 
the increase of the dynamical friction time-scale with host galaxy 
mass, implying that the orbital decay of GCs and their potential 
merging beco mes less efficient in more massive galaxies. Indeed, 



iBekkil ( l2010h identify a disc mass of about lO' Mq as the upper 
limit for efficient GC migration due to dynamical friction and cor- 
responding nGC growth, in good agreement with the data. 

The dwarf galaxy stellar mass range M*.gai = lO'-lO' Mo 
is of particular interest because the stellar halo of the Milky Way 



^ Unless massive clusters are more rapidly disrapted than lo w-mass ones, 
contrary to what is shown by theory and observations (e.g. iLamers et alJ 



l2010l ; lBastian et alj20inih . 



is thoug ht to have formed by the accretion of galaxies with such 
masses jCooper et allbOIOl) . If we account for the disruption his- 
tory of GCs, this can be used to give a first-order estimate for the 
fraction of the current metal-poor Galactic GCs have formed in a 
similar way to YMCs, and which fraction formed as an nGC be- 
fore they were accreted. The masses of outlier nGCs in Fig. |2] are 
about an order of magnitude higher than other GCs in current dwarf 
galaxies, but if they were accreted earlier they may not have had the 
time to grow as much. Taking an intermediate mass difference of a 
factor of five, Eq.|5]shows that the disruption rate of nGCs is about 
three times lower than for GCs. We can then use the difference in 
the specific frequency of GCs between dwarf galaxies and Milky 
Way-mass galaxies in the Virgo Cluster (number of metal-poor 
GC s per unit galaxy stellar mass Tbiuc, with T^^"f / Tj^!^^^ ~ 0.2, 
see IPeng et ^l2008ll to get a rough idea of the total amount of 
disruption during the assembly of the Galactic GC system from ac- 
creted dwarf galaxies. Combining this with the initial nGC fraction 
of 10%, we find that the current fraction of metal-poor Galactic 
GCs that was formed differently than YMCs could well be as high 
as 35%. This fraction is an upper limit because it assumes that 
all current metal-poor Galactic GCs were once the most massive 
cluster of a dwarf galaxy, which need not be the case. Since metal- 
poor GCs constitute about half of the Galactic GC population, the 
Galaxy-wide upper limit to the number of GCs with a nuclear ori- 
gin is 15-20%. 

The fraction of GCs with a nuclear origin is interesting 
in view of the several Galactic GCs with chemical signatures 
that are strongly in consistent with single stellar populations (e.g. 
iGratton et alj 2004]). If these chemi cal anomalies can be traced to 
former nGCs 1 Carretta et alj I2OI0I) . the observed distinct stellar 
populations would require them to be the product of one or two 
GC-GC mergers in the centre of their host galaxies instead of a 
continuous accretion. This is not unlikely, because dynamical fric- 
tion is most efficient for the most massive clusters. According to 
current data, the fraction of Galactic GCs with the light element 
abundance anomalies that indicate multiple stellar populations is 
much larger than 35% iC aloi & D'Antona 201 1) , casting doubt on 
the feasibility of explaining them by nGCs. Instead, it suggests that 
light element abundance anomalies should also be present in the 
most massive YMCs in the nearby universe - depending on how 
long it takes for a second generation of stars to form. However, 
the fraction of GC s with Fe abundance variations is smaller (e.g. 
iGratton et alj2004h . and could possibly be consistent with the esti- 
mated fraction of nGCs. 



We have shown that the majority (65-90%) of the metal-poor 
GCs in dwarf galaxies are consistent with the Schechter-type ICMF 
of YMCs, which is required for their formation mechanisms to have 
been the same. A log-normal ICMF is ruled out with 99.6% sig- 
nificance. However, in about 10% of the dwarf galaxies between 
M*_gai = lO^-lO' M0 the masses of the most massive GCs are 
at least 2a too high to be consistent with YMCs. These outliers 
are generally nGCs, of which the masses have increased due to 
the merging with other clusters and further gas accretion, allowing 
them to exceed the mass-scales set by the star formation process 
(and thus the ICMF). Their elevated masses imply a higher chance 
to survive any tidal disruption during and since their accretion onto 
the Milky Way than 'normal' GCs, suggesting that the current frac- 
tion of metal-poor GCs with a nuclear origin will have risen above 
10%, potentially up to 35%. 



6 /. M. D. Kruijssen and A. P. Cooper 
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